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Abstract 

CN This paper is dedicated to the consistency of systemic risk mea- 

^ sures with respect to stochastic dependence. It compares two alter- 

^ native notions of Conditional Value-at-Risk (CoVaR) available in the 

current literature. These notions are both based on the conditional 
distribution of a random variable Y given a stress event for a random 
variable X, but they use different types of stress events. We derive 
representations of these alternative CoVaR notions in terms of cop- 
ulas, study their general dependence consistency and compare their 
performance in several stochastic models. Our central finding is that 
conditioning on X > \aRa{X) gives a much better response to de- 
pendence between X and Y than conditioning on X = VaRa(X). We 
prove general results that relate the dependence consistency of CoVaR 
using conditioning on X > VaRa(X) to well established results on con- 
cordance ordering of multivariate distributions or their copulas. These 
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Q results also apply to some other systemic risk measures, such as the 

Marginal Expected Shortfall (MES) and the Systemic Impact Index 
(SII) . We provide counterexamples showing that CoVaR based on the 
stress event X = YaRa{X) is not dependence consistent. In particu- 
lar, if {X, Y) is bivariate normal, then CoVaR based on X = YsiRa{X) 
^ is not an increasing function of the correlation parameter. Similar is- 

sues arise in the bivariate t model and in the model with t margins and 
a Gumbcl copula. In all these cases, CoVaR based on X > YsiRa(X) 
is an increasing function of the dependence parameter. 
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1 Introduction 



The present paper studies the notion of Conditional Value-at-Risk (CoVaR) 
introduced by Adrian and Brunnermeier [2] as a dependence adjusted version 
of Value-at-Risk ( VaR) . The general idea behind CoVaR is to use the condi- 
tional distribution of a random variable Y representing a particular financial 
institution (or the entire financial system) given that another institution, 
represented by a random variable X, is in stress. CoVaR represents one of 
the major threads in the current regulatory and scientific discussion of sys- 
temic risk, which significantly intensified after the recent financial crisis. The 
current discussion on systemic risk measurement is far from being concluded, 
and the competing methodologies are still under development. In addition to 
systemic risk measures [cf. 2, 4, 14, 15, 1, 29, 16], related topics include the 
structure of interbank networks, e.g., [8, 11], models explaining how systemic 
risk is created, e.g., [10, 17], and attribution of systemic risk charges within 
a financial system, as discussed in [26, 25]. 

Our contribution addresses the consistency of systemic risk measures with 
respect to the dependence in the underlying stochastic model. In the case 
of CoVaR we give a strong indication for the choice of the stress event for 
the conditioning random variable X. There are two alternative definitions 
of CoVaR in the current literature. The original definition in [2, 3, 4] is 
derived from the conditional distribution of Y given that X = VaRo,(X). 
The second one uses conditioning on X > VaRQ,(X). This modification 
was proposed by Girardi and Ergiin [14] to improve the compatibility of 
CoVaR with non-parametric estimation methods. For similar reasons, such as 
continuity and better compatibility with discrete distributions, conditioning 
on X > VaRa(X) was also favoured by Klyman [20] for both CoVaR and 
Conditional Expected Shortfall (CoES). Finally, it is remarkable that most 
competitors of CoVaR [cf. 15, 1, 29, 16] use conditioning on X > VaRQ,(X) as 
well. This approach goes in fine with the general concept of stress scenarios 
discussed in [6] . 

Our results show that conditioning on X > VaRQ,(X) has great advan- 
tages for dependence modelling. We prove that this modification of CoVaR 
makes it response consistently to dependence parameters in many important 
stochastic models, whereas the original definition of CoVaR fails to do so. 
The counterexamples even include the bivariate Gaussian model, where the 
original CoVaR is decreasing with respect to the correlation p := corr(X, Y) 
for p > 1/v^- Thus, CoVaR based on {X = VaRQ,(X)} fails to detect sys- 
temic risk when it is most pronounced; and we also found this kind of incon- 
sistency in other examples. On the other hand, our findings for the modified 
CoVaR relate its dependence consistency to concordance ordering of multi- 
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variate distributions or related copulas. This may explain the comparative 
results in [13]. where CoVaR^ stood somewhat apart from its competitors. 
Moreover, it gives the the modified notion of CoVaR a solid mathematical 
basis. 

Besides CoVaR, we also discuss extensions to Conditional Expected Short- 
fall (CoES). It turns out that the dependence inconsistency or dependence 
consistency of the alternative CoVaR notions is propagated to the corre- 
sponding definitions of CoES. The dependence consistency results for CoVaR 
and CoES based on the stress scenario X > VaRa(X) also apply to the 
Marginal Expected Shortfall (MES) defined in [1] and to the Systemic Im- 
pact Index (SII) introduced in [29]. 

The paper is organized as follows. Basic notation and alternative defini- 
tions of CoVaR and CoES are given in Section 2. In Section 3 we present the 
general mathematical results, including representations of CoVaR in terms 
of copulas and consistency of the modified CoVaR or CoES with respect to 
dependence characteristics. Section 4 contains a detailed comparison of the 
original and the modified CoVaR in three different models: the bivariate nor- 
mal, the bivariate t distribution, and a bivariate distribution with t margins 
and a Gumbel copula. Conclusions are stated in Section 5. 

2 Basic definitions and properties 

Let X and Y be random variables representing the profits and losses of two 
financial institutions, such as banks. Focusing on risks, let X and Y be 
random loss variables, so that positive values of X and Y represent losses, 
whereas the gains are represented by negative values. 

The issues of contagion and systemic stability raise questions for the joint 
probability distribution of X and Y: 

Fx,Y{x,y) ■.= V{X<x,Y <y). 

The corresponding marginal distributions will be denoted by Fx and Fy- 
Provided a method to quantify the loss or gain of the entire financial system, 
Fx,Y can also represent the joint loss distribution of a bank X and the system 
Y.' 

In the current banking regulation framework (Basel II and the so-called 
Basel 2.5), the calculation of risk capital is based on measuring risk of each 
institution separately, with Value-at-Risk (VaR) as a risk measure. The 
Value-at-Risk of a random loss X at the confidence level a e (0, 1) is the 
a-quantile of the loss distribution Fx [cf. 21, Definition 2.10]. That is, 

VaR„(X) = F^(a) 



3 



where F^iy) '■— ini{x e R : Fx{x) > y} is the generalized inverse of Fx- 
The most common values of a are 0.95 and 0.99. 

For continuous and strictly increasing Fx the generalized inverse F^ co- 
incides with the inverse function of Fx- In this case one has VaRa(X) = 
F^^{a) for a G (0, 1). For a thorough discussion of generalized inverse func- 
tions we refer to [12]. 

In the present paper we discuss two alternative approaches to adjust 
VaR to dependence between X and Y. This is achieved by conditioning the 
distribution of F on a stress scenario for X. These two notions appear in 
the recent literature under the name Conditional Value-at-Risk (Co VaR), 
but they use different kinds of stress scenarios. The original notion was 
introduced by Adrian and Brunnermeier [2, 3, 4] and will henceforth be 
denoted by Co VaR". The alternative definition was proposed by Girardi 
and Ergiin [14]. We denote it by Co VaR. 

Definition 2.1. 

CoVaR=^(r|X) := VaR;3(>"|^ = VaR„(X)); 
CoVaRa,/3(F|X) := VaR^(F|X > VaR«(X)). 

The computation of CoVaR" requires the knowledge of FY\x=YaRa{x)- If 
Fx,Y has a density /x,y, then fx{x) = J^^ fx,Y{x,y)dy is a density of Fx, 
and 

, , X!^/x,y(VaR«(X),i)di 
Fy|x=VaR„(x)(y) j^(VaR,(;f)) ' 

provided that /x(VaRa(X)) > 0. In some models, such as elliptical distri- 
butions, FY\x=YaRaix) is known explicitly. In general, however, computation 

of i^y|x=VaRa(x) requires numerical integration. 

Conditioning on X > VaRQ,(X) is less technical. The definition of VaRQ,(X) 
implies that P(X > VaRQ(X)) > 1 — a, so that elementary conditional prob- 
abilities are well defined. In particular, if Fx is continuous, then 

P(y <Z/,X> VaR.(X)) 

^Y\X>Y3.RociX)[y) — : ■ 

1 — a 

Moreover, conditioning on events with positive probabilities is advantageous 
in statistical applications, including model fitting and backtesting. This is 
the major reason why the original notion of CoVaR^ was modified to CoVaR 
in [14]. 

A straightforward extension from CoVaR to Conditional Expected Short- 
fall (CoES) is based on the representation ES^(y) = J^VeiRt{Y)dt. 
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Definition 2.2. 

1 

CoESaAY\X) := — CoV8.R^,t(Y\X)dt, (1) 

CoES=^(y|X) := ^ CoVaR=„,*(y|X)c/t (2) 

Remark 2.3. (a) In precise mathematical terms, CoVaR^^ and CoVaRa,/? 
are the /3-quantiles of the conditional distributions -FV|x=VaRa(x) and 

FY\X>YaRa(xy- 

CoVaR=^(r|X) = F^x=vaR„(x)(/3); 
CoVaR„,^(y|X) = F^x>vaR„(x)(^)- 

(b) In [2, 4, 18, 14], the authors work with a common confidence level for 
X and F, i.e., in the special case a = (3. Similarly to the notation used 
there, we will omit /3 ii /3 — a and write CoVaRa instead of CoVaRa^a 
if it does not lead to confusion. However, the definition of CoES needs 
separate confidence levels for X and Y in the integrand CoVaRa,t{Y\X). 

(c) Since CoESQ,,;3(y|X) = ESp{Z) for a random variable Z ~ FY\x>YaRc{x), 
the coherence of ES in the sense of [5] is inherited by CoESq,^ for all 
a, /3 e (0, 1). The central point here is subadditivity, which is understood 
as 

CoES«,^(r + Y'\X) < CoESaAY\X) + CoES„,^(F'|X) 

for any random variables (F, Y', X) defined on the same probability 
space. 

(d) In [2, 4], CoES is defined as E[Y\Y > CoVaR=„(r |X)]. Note that 
this definition replaces the stress scenario {X = VaRa(X)} by {Y > 
CoVaR^(^(F|X)}, which is not related to X directly. Compared to 
CoES^_^(y|X), this definition is quite unnatural. Moreover, it does not 
guarantee coherence, which is the central property of Expected Shortfall. 

(e) The notion of Marginal Expected Shortfall (MES) introduced in [1] is 
closely related to CoVaR and CoES. It is defined as 

MESa{Y\X) := E[Y\X > VaR^{X)] 

where X := Yli=i is the financial system and Y ■.= Yi for some i is an 
institution. The idea behind MES is to quantify the insurance premia 
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corresponding to bail-outs which become necessary when the entire finan- 
cial system is close to a collapse. The major economic difference between 
MES and CoVaR is the role of X and Y. With MES, the conditioning 
random variable X is the system, and the target random variable Y is 
a part of the system. In the original work on CoVaR, Y is the system, 
and X is a part of it. 

On the mathematical level, MES and CoVaR or CoES are quite close to 
each other. It is easy to see that 

MESaiY\X)= /'F^x>vaR„(x)Wrf^= r CoYaR^,tiY\X)dt. 
Jo Jo 

In view of (1), one could also write MESa{Y\X) = CoESa,o(^|^). 

(f) In [20], CoVaRo,,/? and CoESa,/? in the sense of Definitions 2.1 and 2.2 
are called DistVaR and DistES. Besides the different naming, the defi- 
nitions are essentially the same, and these notions are also compared to 
CoVaR^ and CoES". However, the comparison in [20] is concentrated on 
general representations, compatibility with discrete, e.g., empirical, dis- 
tributions, and the behaviour in the bivariate Black-Scholes model. As 
far as we are aware, a study of consistency with respect to dependence 
parameters has been missing so far. 

The introduction of CoVaR" in [2] aims not at CoVaR" itself, but at 
the contribution of a particular financial institution to the systemic risk. 
In [2] , CoVaR" is used to construct a risk contribution measure that should 
quantify how a stress situation for an institution X affects the system (or 

another institution) Y. In [2], the authors propose — 1 as a 

systemic risk indicator. In [3], the systemic risk measure is modified to 

ACoVaR=^(F|X) := CoVaR" ^(F) - VaR^iY). (3) 

In [4], the centring term VaR^(y) representing the risk of Y in an unstressed 
state is replaced by the conditional VaR of Y given that X is equal to its 
median: 

A'"'''^CoVaR=^(r |X) := CoVaR=^(y|X) - VaR^(r|X - med(X)) (4) 

to remedy some inconsistencies observed in a comparison of CoVaR" across 
different models. 

Unfortunately, the centring in (3) is not the only reason why ACoVaR" 
can give a biased view of dependence between X and Y. The results presented 
below demonstrate that there is a more fundamental issue that cannot be 
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solved by modifying ACoVaR= to A"^^^CoVaR= or taking any other centring 
term. The primary deficiency of ACoVaR^ is that the underlying stress 
scenario X = VaRa(X) is too selective and over-optimistic. If, for instance, 
Fx is continuous, then P{X = VaRa(X)) = 0, so that this particular event 
actually never occurs. Generally speaking, the ability of CoVaR", ACoVaR", 
or A™^^CoVaR= to describe the influence of X on y strongly depends on 
how well FY\x=V3.Ra{x) approximates Fy\x=x ^ot x > VaRa(X). As shown in 
Section 4, this approximation fails even in very basic models, and it typically 
underestimates the contagion from X to Y. 

3 General results 

We begin with representations of CoVaR" and CoVaR in terms of copulas. 
It is well known that any bivariate distribution function Fx,y admits the 
decomposition 

FxA^,y)=C{Fxix),Fyiy)) (5) 

where C is a probability distribution function on (0, 1)^ with uniform margins 
(cf. [24, 19]). That is, there exist random variables U,V ^ unif(0, 1) such 
that C{u,v) = P{U < u,V < v). The function C is called a copula of 
FxY- If both Fx and Fy are continuous, then C is uniquely determined by 
c(u,v)^Fx,y{F^{u),F^{v)). 

The decomposition (5) yields the following representation of CoVaR~ and 
CoVaR. 

Theorem 3.1. Let {U,V) ~ C where C is a copula of Fx,y- If Fx is 
continuous, then 

(a) CoVaR=^(r|X) = 

(b) CoYaR^Ay\X) = {F^\u>ai^))' Fv\u>a{v) = ^^=g^. 
Proof. Part (a). It is well known that {Fy iU) F^ iV)) ~ Fx,y-i and hence 

FY\x=v.K^ix){y) = P(i^^(^) < y\F^{U) = F-(a)). 

The functions Fy and Fy are non-decreasing and satisfy v < Fy{Fy{v)) and 
F^{FY{y)) < y for aU t; e (0, 1) and y e R. This implies that F^{V) < y 
is equivalent to V < Fyiy). Moreover, continuity of Fx imphes that Fx is 
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strictly increasing, so that F^{U) — F^{a) is equivalent to U — a. This 
yields 

FY\x=Y.R„ix){y) = P{V < FY{y)\U = a) = Fv\u=aiFY{y)), 

and the result follows from the chain rule for the generalized inverse. 
Part(b). Analogously to Part (a), one obtains that 

FY\x>v.R^(x){y) = P{V < FY{y)\U > a) = Fviu>a{FY{y)), 

and hence CoVaR„,^(r |X) = F^{Fv\u>a{/3)). Since {U,V) - C and the 
margins of C are uniform, we obtain that 

,, P{V<v,U>a) v-C{a,v) 

Fv\U>a[V) = 



P(U>a) I- a 

□ 

Theorem 3.1(a) provides a link between the ordering of CoVaR and the 
notion of concordance ordering. 

Definition 3.2. [cf. 22, Definition 3.8.1] Let (X, Y) and {X\ Y') he bivariate 
random vectors with Fx — Fx' and Fy — Fy'. Then (X, Y) is smaller than 
{X' ,Y') in concordance order ({X,Y) -< {X' ,Y') or, equivalently, Fx,y di 

Fx',Y') if 

yx,yeR F{X <x,Y <y) <F{X' <x,Y' <y)- 

Remark 3.3. The following equivalent characterizations of {X,Y) ^ [X',Y') 
will be used in in the sequel: 

(a) P(X >x,Y>y)< P{X' > xX > y) for all x, y e R; 

(b) C d C for the copulas of Fx,y and Fx',y' if the margins are continuous; 

(c) Ef{X,Y) < Ef{X',Y') for all supermodular functions / : ^ M, i.e., 
for all / satisfying 

fix + e,y + S) + fix, y) > fix + £, y) + /(x, y + 5) 

for all x, y e R and all £, 5 > 0. This order relation is called supermodular 
ordering (:<sm)- 

For proofs and further alternative characterizations we refer to [22, Theorem 
3.8.2]. 



The central theoretical result of the present paper is the following. 

Theorem 3.4. Let {X, Y) and {X', Y') be bivariate random vectors with 
copulas C and C , respectively, and assume that Fy = Fy- 

(a) If Fx and Fx' are continuous, then C <C' implies 

Va, ^ e (0, 1) C6VbR^,p{Y\X) < CoVaRc,,^(y'|X'). (6) 

(b) If Fx, Fx', Fy, and Fyi are continuous, then (6) implies C ^ C . 

Remark 3.5. Note that Theorem 3.4 does not need Fx — Fx'- The only 
assumption on the conditioning random variables X and X' is that they are 
continuously distributed. 

Proof Part (a). Let {U,V) ~ C and {U' ,V') ~ C . As F^ and F^ are 
non-decreasing, Theorem 3.1(b) reduces the problem to 

Va,/3 G (0, 1) F^f;>J/3) < (7) 

Moreover, it is well known that for any distribution functions G and H the 
ordering G^{y) < H*"{y) for all y G (0, 1) is equivalent to G{x) > H{x) for 
all a; G M. Thus it suffices to show that 

ya,v G (0, 1) Fv\u>a{v) > Fv'\U'>a{v)- 

The representation of Fy\ua{v) in Theorem 3.1(b) reduces this to C{a,v) < 
C'{a,v) ior all a,v, which is precisely C :< C'. 

Part (b). Combining (6) with Theorem 3.1, one obtains 

\faj F-(F^^>J/5))<Ff(F-|^,>J/3)). (8) 

As Fy is continuous, F-^ is strictly increasing. Therefore (8) implies (7), 
which is equivalent to C ^ C". □ 

Theorem 3.4 can be applied to various stochastic models. We start with 
elliptical distributions. This model class includes such important examples 
as the multivariate Gaussian and the multivariate t distributions. Since 
CoVaRQ,^^(X|F) considers two random variables and multivariate elliptic- 
ity implies bivariate ellipticity for all bivariate sub- vectors, we restrict the 
consideration to the bivariate case. 

A bivariate random vector {Xi,X2) is elliptically distributed if 

{X,Yy = + RAW^ 
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where = {lix^yt-v) e and A e R^^^ are constant, W = (l^i,W^2) is 
uniformly distributed on the Euchdean unit sphere {x G : ||a;||2 = 1}, and 
i? is a non-negative random variable independent of W^. If ER < oo, then 
= EX and = EY. The ellipticity matrix S := A is unique except 
for a multiplicative factor. The covariance matrix of (X, Y) is defined if and 
only if ER^ < oo, and this matrix is always equal to cE for some constant 
c > 0. Thus, rescaling R and A, one can always achieve that 



where, if defined, ax = var(X), cry = var(F), and p = coyy{X,Y). In the 
following we will always assume this standardization of E and denote the 
bivariate elliptical distribution with location parameter p — [pxiPy] and 
ellipticity matrix S by E{fi,Ti, R). 

If {X,Y) ~ SIpjUjR) with continuous marginal distributions, then the 
copula C of {X, Y) is uniquely determined. Copulas of this type are called 
elliptical copulas. The invariance of copulas under increasing marginal trans- 
forms implies that C depends only on the parameter p of E and on the 
distribution of R. Thus p is the natural dependence parameter for a bivari- 
ate elliptical copula C, whereas the distribution of R specifies the type of 
the copula, such as Gaussian or t. We will call elliptical copulas C and C 
of same type if the corresponding elliptical distributions have identical radial 
parts R^ R'. 

The following theorem states monotonicity of CoVaR with respect to 
the dependence parameter p if (X, Y) is elliptically distributed or has an 
elliptical copula. In particular, it applies to bivariate Gaussian or bivariate t 
distributions, and also to bivariate distributions with Gaussian or t copulas. 

Theorem 3.6. (a) Let (X, Y) ~ £{p, E, R) and {X', Y') ~ £{p', S', R) with 
continuous Fx and Fx'. If Py < Pv CLnd uy = cry, then p < p' implies 



(b) Let (X,F) ~ S{p,^,R) and {X',Y') ~ S{p',J:',R) with continuous Fx 
and Fx'. If py < py and ay < ay, then p < p' implies 

Va e (0, 1) e [^0, 1) CoVaRc,,^(y|X) < CoVaR«,^(y'|X') 
with f3o := ^'^^"^"'"'"^^^ where C is the copula of (X, F). 

(c) Let Fxy and Fx'.y have elliptical copulas of same type with dependence 
parameters p and p' , respectively. If Fx and Fx' are continuous and 
Fviy) > FY'{y) for all |/ G M, then p < p' implies (6). 





(6). 
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Remark 3.7. (a) The assumption Fy > Fy obviously includes the case of 
identical margins Fy = Fy' , which is the natural setting for studying the 
response of CoVaR to dependence parameters. 

(b) It is easy to see that the lower bound /3o in Theorem 3.6(b) is decreasing 
in p. In particular, one has Pq < 1/2 for p > 0. This guarantees that 
CoVaR„,;3(F|X) < CoVaR„,/3(r'|X') for G [1/2,1), which is fully 
sufficient for assessing dependence between rare events. 

Proof of Theorem 3.6. Part (a). It is obvious that CoVaKa.B{c + ^|-^) = 
c + CoVaRc,^(F|X). Hence, as fiy < pyi, it suffices to consider /ly = /lyi, 
so that we have Fy — Fyi . Since the case cry = is trivial, we only need to 
consider cry > 0. 

The continuity of Fx yields ax > 0, and as {X,Y) is elliptically dis- 
tributed, we have Y = —X. Hence Fy is continuous as well, and therefore 
the copulas C and C of (X, Y) and (X , Y ) are uniquely defined. 

According to Theorem 3.4(a), it suffices to show that p < p' implies 
C ^ C. This is equivalent to £ (0, 0, r(p), i?) ^ £"(0, 0, r(p'), 7?) for p < p' 

and r(p) = ^ ^ 1 '^^^^ ordering result is proven in [9]. In the bivariate 

Gaussian case it is also known as Slepian's inequality [cf. 27, Theorem 5.1.7]. 

Part (b). Without loss of generality we can assume that py = py and 
cry > 0. Part (a) gives us C :< C" and hence (7). Moreover, cry < cry/ implies 
that F^{t) < F^{t) for t e [1/2,1). Hence, according to Theorem 3.1(b), 
it suffices to verify that Fyp^^{P) > 1/2. This inequality is equivalent to 

P > Fviu>a{l/2) = (3o. 

Part (c). According to Part (a), we have C ^ C and hence (7). Since 
Fy{y) > Fy\y) for all |/ G M is equivalent to Fy{y)^{t) < F^jt), Theo- 
rem 3.1(b) yields 

CoVaR„,^(F|X) < F^,{F^^u>am < CoVaR„,^(r|X'). 

□ 

A very popular copula model is the Gumbel copula. In the bivariate case 
it is defined as 

CeM = exp (- {{-loguf + {-logvff") . (10) 

The dependence parameter 9 assumes values in [1, oo], whereas 9 — 1 and 9 — 
oo refer to Ci{u,v) := uv (independence copula) and Coo{u,v) := mm{u,v) 
(comonotonicity copula). As shown in [28], 9 < 6' implies Cq Cgi and 
hence Cg ^ Cgi (cf. Remark 3.3(c)). This immediately yields the following 
analogue of Theorem 3.6(c). 
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Corollary 3.8. Let {X,Y) and {X',Y') have Gumbel copulas with depen- 
dence parameters 9 and 6' , respectively. If Fx and Fx' are continuous and 
Fyiv) > Fy'iy) for all yeR, then 9 < 9' implies (6). 

Remark 3.9. Corollary 3.8 also holds for Galambos copulas with dependence 
parameters 9 < 9'; see [28] for Cq Cqi in this case. 

The monotonicity of CoESa,i3{X,Y) with respect to dependence param- 
eters follows from the integral representation (1). 

CoroIIeiry 3.10. Suppose thatE\Y\ andE\Y'\ are finite. 

(a) If{X,Y) and {X',Y') satisfy the assumptions of Theorem 3.6(a) or (c), 
or those of Corollary 3.8, then 

V«,/3 e (0, 1) CoES„,^(r|X) < CoES«,^(F'|X'). (11) 

(b) If {X,Y) and {X',Y') satisfy the assumptions of Theorem 3.6(h), then 

e (0, 1) e [^0, 1) CoESa,p{Y\X) < CoES«,^(y'|X'). 

wuh = 

We conclude this section by relating the results obtained here to another 
systemic risk measure. 

Remark 3.11. (a) Corollary 3.10(a) also applies to the Marginal Expected 
Shortfall from [1]. Setting /3 = in (11) and applying Remark 2.3(e), 
one obtains MES«(F|X) < MES„(F'|X') for all a. 

(b) In [29] , the Systemic Impact Index (SII) of an institution Yi is defined as 



d 



Yi > VaR«(Fi) 



Slli(a) := E l{Yj > VaR^{Yj)} 

\j=i 

= 1 + ^ F{Yj > Va.Ra{Yj)\Yi > VaR«(FO)- 
It is easy to see that (6) is equivalent to 

P(y > VaR^(y)|X > YaiRp(X)) < P(Y' > YaRf)(Y')\X' > VaR^(X')) 

for all a, (3. Thus, for Y = Yj and X = Yj, the assumptions of The- 
orems 3.4(a) and 3.6 also imply dependence consistency of the single 
conditional default probabilities P{Yj > YaRa{Yj)\Yi > YaRa{Yi)). 
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4 Examples 



In this section we compare CoVaR and CoVaR^ in three different models: 
the bivariate Gaussian, the bivariate t, and the bivariate distribution with a 
Gumbel copula and t margins. 



4.1 The bivariate Gaussian distribution 

It is well Icnown that the bivariate Gaussian distribution is elliptical. Hence 
Theorem 3.6(a) guarantees that CoVaR is an increasing function of the cor- 
relation parameter p. Moreover, CoVaR^ can be calculated explicitly in this 
case, so that it is particularly easy to compare CoVaR with CoVaR". 



Computation of CoVaR" 

Let (X, Y) ~ A/'(yU, S) with mean vector p, = {px, Py) and covariance matrix 
S as in (9). As for all bivariate elliptical models, the dependence between 
X and Y is fully described by the correlation parameter p. An appealing 
property of the bivariate normal distribution is the interpretation as a linear 
model. Indeed, {X, Y) ~ Afip, S) is equivalent to 

I^ = ,^^^ + Vr37z, (12) 

cry ax 

where X ~ M{px^ (^"x) and Z ~ A/'(0, 1), independent of X. 

Due to X ~ J\f{px,o'x) we have VaRa(X) = px + crx$~^(c(;), where $ 
is the distribution function of A/'(0, 1). Substituting X — YaRa{X) in (12), 
one obtains 

Y^Py + (Ty (p$"'(Qi) + Vl - P'^) ■ 

This shows that JC{Y\X = YaRa{X)) = Af{p, a'^) with fx ^ py + aYP^~^{a) 
and a — cry^/l — Hence we obtain that 

CoVaR= .(F|A:) = VaR^(F|X = VaR„(X)) =p + a(!>-\(3) 



Py + CJY [p^-\a) + $-i(/3)Vl-p') ■ (13) 



Computation of CoVaR 



To compute CoVaR, we use the copula representation from Theorem 3.1(b). 
From Y ~ J\f{pY.a^) one obtains that Fy^{v) = py + cry^^^iv) for v G 
(0,1). Moreover, the copula of {X,Y) ~ J\f{p,T,) is the Gauss copula Cp 
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with dependence parameter p. For p = it is the independence copula, 
Cq{u,v) — uv, and for p 7^ it has the following representation: 



C,{u,v) = Fx,y{Fx{u),F^{v)) 



J -00 J -00 27rv/l - p2 ^ V 



2^ 



ds2dsi. 

(14) 



27rVW' V 2(1 -p 

Applying Theorem 3.1(b), we obtain 

CoVaR«,^(y) = py + aY^-\F-^^^J/3)) 

where Fv\u>a{v) = ^^—jz^^- The values of CoVaR can be obtained by nu- 
merical integration of (14) and numerical inversion of the function Fv\u>a{i^)- 
An alternative method to compute CoVaR is the numerical computation 
and inversion of the function 

FY\x>YaRa(x){t) = I I fx,Y{x,y)dxdy, (15) 

^ ~ J-00 JVaRa(X) 

where fx,Y is the joint density of X and Y. Depending on the application, 
each method has its advantages. Whilst (15) is more direct and hence faster 
for numerically tractable /x,y, the conditional copula values obtained in (14) 
can be re-used with other marginal distributions. 

Monotonicity in p 

As bivariate Gaussian distributions are elliptical. Theorem 3.6(a) guarantees 
that CoVaR is always increasing in p. However, this is not the case for 
CoVaR". Partial differentiation of (13) in p yields 

9,CoVaR=,(y|X) = ay (^^-\a) - , (16) 

which is positive if ^~^{a)^/Y^^^ > p^~^(^) and negative if ^~^{a)^/Y^^^ < 
p$~^(/3). Besides the degenerate case a = /3 = 1/2 with constant CoVaR^^, 
there are 4 cases depending on the signs of ^~^{a) and $~^(/3): 

(i) If a > 1/2 and /3 > 1/2, then CoVaR=^(y|X) is increasing in p for 
p < po := , ^"^^ and decreasing for p > pn. 

(ii) If a > 1/2 and (3 < 1/2, then CoVaR^^^(F|X) is increasing in p for 
p > — Po and decreasing for p < — po. 
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CoVaR as a function of p CoVaR as a function of p 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 

p p 



Figure 1: CoVaR=(F|X) and CoVaR«(y|X) (i.e., with /3 = a) in the bivari- 
ate normal model as functions of p. 

(iii) If a < 1/2 and (3 > 1/2, then CoVaR^^^(F|X) is increasing in p for 
p < — po and decreasing for p > — po. 

(iv) If a < 1/2 and /3 < 1/2, then CoVaR^^^(y |X) is increasing in p for 
p > Po and decreasing for p < po. 

Thus CoVaR^ is monotonic with respect to p only in degenerate cases. In 
particular, in the most important case a, /3 G (1/2, 1), CoVaR" is decreasing 
for p > Po, which means that CoVaR" fails to detect dependence where it 
is most pronounced. In the special case a = /3, the critical threshold po is 
always equal to 1/ \/2. 

A graphic illustration to this fact is given in Figure 1, showing CoVaR^(y |X) 
and CoVaRQ,(F|X) for p > 0.2 and a = (3 assuming values 0.90, 0.95, or 0.99. 
The short writing CoVaR^ refers to CoVaR^^; analogously, CoVaR^ denotes 
CoVaRfj Q,- This notation was used in the original definitions of CoVaR" and 
CoVaR, which were restricted io a = (5 (cf. Remark 2.3(b)). For the sake of 
simplicity we set /xy = and ay = 1. These parameters have no influence 
on the decreasing or increasing behaviour of CoVaR or CoVaR" as functions 
of p. 

Normalized values of CoVaR and CoVaR" 

The relative impact of a stress event for X on the institution Y can be quanti- 
fied by the ratio CoVaR" ^ (F | X) / VaR„(F) or by CoVaR„,/3(r |X)/ VaR«(y). 
A similar of systemic risk indicator was proposed in [2] . Figure 2 shows these 
ratios for a = P and p = as functions of a. The different line types in the 
plots correspond to p = 0.5, 0.7, and 0.92. The ratios CoVaR^(y|X)/VaRQ,(F) 
are constant, which is also easy to see from (13). The interesting part here 
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Ratio CoVaR°A/aR as a function of a and p Ratio CoVaRA/aR as a function of a and |3 

2.0 -| p 

0.5 

0.7 

1.8 0.9 

1.6 - 



1.4 - 



1.2 - 

T 1 1 1 1 r 

0.90 0.92 0.94 0.96 0.98 1.00 0.90 0.92 0.94 0.96 0.98 1.00 



Figure 2: Bivariate normal model with /iy = 0: Ordering of the ratios 
CoVaR=(F|X)/ VaR„(F) and CoVaR„(F|X)/ VaR«(F) for different a. 

is the ordering of the lines for different p. In case of CoVaR^, the line for 
p = 0.7 is above the two others, illustrating that the inconsistency issue is 
common to all a G (1/2,1). The plot of CoVaR„(F|X)/ VaRc,(F) shows 
correct ordering for all a, as guaranteed by Theorem 3.6(a). Another obser- 
vation one can make here is that CoVaIia{Y\X) / VaRQ,(y) is decreasing in a. 
This, however, is a model property that seems to be related to the light tail 
of the normal distribution. In heavy-tailed models considered in Sections 4.2 
and 4.3 the ratio CoVaRQ,(y|X)/ VaRQ,(F) is increasing in a. 

Backtesting and violation rates 

The results above show that CoVaR reflects the dependence between X and 
Y much more consistently than CoVaR^. An intuitive and very general 
explanation to this fact is that conditioning on X > VaRa(X) corresponds to 
a reasonable "what if question, whereas conditioning on X = VaIia{X) does 
not. Indeed, the scenario {X > VaRQ(X)} includes all possible outcomes for 
X if X is stressed, whereas the scenario {X = VaRa(X)} selects only the 
most benign of them. 

In backtesting of VaR one expects that X exceeds VaRQ(X) with prob- 
ability not larger than 1 — a. Abbreviating "Conditional VaR", the term 
CoVaR^^^ suggests that Y exceeds CoVaR^^^(F|X) with conditional prob- 
ability 1 — /3 or less, given that X is stressed. The definition of Co VaR 
understands stress of X as {X > VaRa(X)}, so that the expected violation 
rate for CoVaRa,i3 under this stress scenario is equal to 1 — /3. In contrast to 
that, CoVaR^ is designed to have the violation rate 1 — /3 under the less nat- 
ural and more optimistic scenario {X = VaRQ(X)}. As a consequence, the 
violation rates for CoVaR^o backtesting experiments based on the natural 
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stress scenario {X > yaRa{X)} are significantly higher than 1 — /3. 

This issue is illustrated in Table 1. The underlying Monte Carlo exper- 
iment generates an i.i.d. sample (Xj.Kj) ~ A/'(0,S) for i = l,...,n and 
counts the joint exceedances {Yi > CoVaR^^(y |X), > YaRa{X)}. The 
CoVaR" violation rate for the stress scenario {X > VaRa(X)} is the ratio of 
the joint excess count and the count of the excesses {Xi > VaRQ,(X)}. The 
violation rate for CoVaR is obtained analogously from the number of joint 
exceedances {Yi > CoVaR„jY\X),Xi > VaRa(X)}. We chose n = 10^ and 
a,/3 being either 0.95 or 0.99. 

It is remarkable that the violation rate for CoVaR^ increases with p. 
This demonstrates that the underestimation of risk by CoVaR^ is most pro- 
nounced in case of strong dependence and, hence, high systemic risk. 



Bound p = p = 0.2 p = 0.5 p = 0.7 p = 0.9 



CoVaRo=95,o.95(>^l^) 


0.0503 


0.0601 


0.0857 


0.1229 


0.2520 


CoVaRo=99,o.99(r|X) 


0.0099 


0.0124 


0.0189 


0.0292 


0.0875 


CoVaRo=95,o,99(>^l^) 


0.0101 


0.0130 


0.0213 


0.0375 


0.1224 


CoVaRo=99,o.95(^l^) 


0.0500 


0.0588 


0.0785 


0.1045 


0.2053 


CoVaRo.95,0.95(l^|^) 


0.0503 


0.0500 


0.0503 


0.0495 


0.0499 


CoVaRo.99,0.99(>^|X) 


0.0099 


0.0101 


0.0104 


0.0099 


0.0098 


CoVaRo.95,0.99(>^|X) 


0.0101 


0.0102 


0.0102 


0.0099 


0.0098 


CoVaRo.99,0.95(>^|^) 


0.0500 


0.0507 


0.0509 


0.0501 


0.0491 



Table 1: Violation rates in the bivariatc normal case. Monte Carlo backtest- 
ing with n = 10^ and a, /3 e {0.95, 0.99} 

A graphical illustration to this issue is given in Figure 3 by bivariate 
normal samples from the simulation study described above. The horizontal 
lines mark the levels of CoVaR=(y|X) and CoVaRa(F|X), and VaR^{Y). 
The vertical lines mark VaRQ,(X). The joint excess counts are the num- 
bers of points above the corresponding horizontal line and on the right 
hand side from the vertical hue marking VaRQ,(X). The sample size is 
n — 2000, which suffices to demonstrate how correlation changes the shape 
of the sample cloud and thus increases the number of the joint excesses 
{Yi > CoVaR=(F|X),X, > VaR«(X)}. 

Risk contribution measures ACoVaR= and A™'^'^CoVaR" 

As mentioned in Section 2, [2] aims not at CoVaR" itself, but at the difference 
between CoVaR" and some characteristic of an unstressed state. The two 
most common definitions of such a risk contribution measure are ACoVaR" 
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Bivariate Normal for p = Q 




Figure 3: Bivariate normal samples (size n = 2000) and the joint excess 
regions in the backtesting experiment for a = /3 = 0.95. 



ACoVaR as a function of p 



A"''"'CoVaR as a function of p 
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Figure 4: ACoVaR„ and A™*^'^CoVaR^ as functions of p in the bivariate 
normal model. 
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and A™^^CoVaR~ (see (3) and (4)). In the bivariate normal case one has 
VaR^(y) = /xy + ay^'^iP), so that (13) yields 



ACoVaR=^(F) = ay (^"^ {a)p + <!>-\/3) (^/T^^. i)) . 

For a = /3 this simplifies to ACoVaR=(F) = crY<^-\a) (^p + y^l - - l). 

Regardless of a and /3, ACoVaR^ inherits the non-monotonicity in p from 
CoVaR^. An illustration to this issue is given in Figure 4, which shows plots 
of ACoVaR= and A™^'^CoVaR= as functions of p for a = ^. 

At a first glance, A™'''^CoVaR" seems to be an improvement because it is 
increasing in p. In fact, A™'^'^CoVaR^ is even linear here. Due to med(X) = 
px, (12) yields F^;,=^ed{x) (Z^) = I^y + ay^l-p^^-\l3). Applying (13), 
one obtains that 

A-^CoVaR=^(y) 

^py + ay ($-'(a)p + $-'(/3) Vl - p") - (py + ay^-\P)^Jl-p'^ 
= (7Y$-^(a)p. (17) 

Thus, in the bivariate normal model, A™'^'^CoVaR^^(y |X) is linear with 
positive slope that depends on p and a, but not on /3. In view of the linear 
structure (12) of the bivariate Gaussian model, this even appears reasonable. 
However, examples in Sections 4.2 and 4.3 show that A™^'^CoVaR^ is not a 
monotonic function of dependence parameters in other models. Thus the ap- 
plicability of A™'^'^CoVaR" is restricted to linear models of type (12), where 
it is superfluous because it carries quite the same information as the correla- 
tion parameter p or the linear regression parameter from the classical Capital 
Asset Pricing Model (the so-called CAPM-^), which is equal to pay /ax in 
the present setting. 

Extension from CoVaR to CoES 

Due to Corollary 3.10(a) we already know that CoESa,/? is increasing in p 
for all a and /3. The special case a = /3 is illustrated in Figure 5, which 
shows that CoES" is not increasing in p. Due to the light tail of the nor- 
mal distribution, these plots are similar to those of CoVaR and CoVaR" in 
Figure 1. A closer look at (2) confirms that the non-monotonicity of CoES" 
in p is inherited from CoVaR". Thus the best possible extension to Condi- 
tional Expected Shortfall based on CoVaR" still fails to reflect dependence 
properly. 
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CoES as function of p 



CoES as function of p 




Figure 5: CoES^(Y\X) and CoESQ(y|X) in the bivariate normal model as 
functions of p. 



4.2 Bivariate t distribution 

The next example we consider is the bivariate t distribution, which is ellip- 
tical, but heavy-tailed. The comparison follows the same scheme as in the 
previous section. A bivariate t distributed random vector with u > degrees 
of freedom (bivariate can be obtained as follows: 

{X,Y) :=(^^,^^) + yZ(^x,f), 

where (X,?) ~ A/'(0, S) and W ~ x'^i^), independent of {X,Y). The 
parameters Px.Py G specify the location of (X, F). For simplicity, we 
consider a centred model with px = /xy = 0. 

It is well known that the bivariate t distribution is elliptical with ellipticity 
matrix S. The corresponding sample clouds have an elliptical shape (cf. 
Figure 8). The second moments of X and Y are finite for i/ > 2, and in 
this case the correlation between X and Y is equal to p. The role of p is 
the same as for all elliptical models: larger values of p increase association 
between large values of X and Y . Analytic expressions for CoVaR^ or CoVaR 
are not feasible in this model, so that computations have to be carried out 
numerically. 



Monotonicity in p 

The behaviour of CoVaR^ and CoVaR as functions of the correlation pa- 
rameter p is shown in Figure 6. Similarly to the Gaussian case, CoVaR is 
increasing in p due to Theorem 3.6(a), whereas CoVaR^ is not. Moreover, 
the relative distance between CoVaR" and CoVaR (as it could be quantified 
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CoVaR as a function of p CoVaR as a function of p 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 



P P 

Figure 6: Bivariate t{3) distribution: CoVaR^ and CoVaRo, as functions of 
the correlation parameter p. 

Ratio CoVaR~A/aR as a function of a and p Ratio CoVaRA/aR as a function of a and p 




n 1 \ \ \ r n \ 1 \ 1 r 

0.90 0.92 0.94 0.96 0.98 1.00 0.90 0.92 0.94 0.96 0.98 1.00 

a a 

Figure 7: Bivariate t{3) distribution with /iy = 0: Ordering of the ratios 
CoVaR=(F|X)/ VaR„(F) and CoVaR„(F|X)/ VaR„(y) for different a. 

by the ratio CoVaR / CoVaR^) is larger than in the Gaussian case. A possible 
explanation to this effect could be the heavy tail of the t{3) distribution. 

Normalized values of CoVaR and CoVaR" 

Figure 7 shows the ratios CoVaR=(F|X)/ VaR^(r) and CoVaR„(F|X)/ VaR^(X) 
as functions of a for selected values of p. This comparison is analogous to 
Figure 2 in the Gaussian case. Similarly to the Gaussian case, the ordering of 
CoVaR^ / VaRa with respect to the dependence parameter p or 6^ is inconsis- 
tent, whereas the ratios CoVaR^/VaRa are ordered correctly for all a: the 
line for the largest p is entirely above the line for the second largest p, etc. 
In contrast to the Gaussian case, these ratios are increasing in a. This could 
be explained by the heavy tail of the t(3) distribution or by the positive tail 
dependence in the bivariate t model. 
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p = p = 0.2 p = 0.5 p = 0.7 p = 0.9 

CoVaR^95^o.95(^l^) 0.1017 0.1213 0.1659 0.2202 0.3638 

CoVaR^gg'o'ggiFlX) 0.0358 0.0433 0.0643 0.0939 0.1909 

CoVaR^g5'o^gg(F|X) 0.0341 0.0429 0.0640 0.0944 0.1954 

CoVaR=gg^o'g5(r|X) 0.1036 0.1229 0.1658 0.2184 0.3546 

CoVaRo.95,o.95(>"|^) 0.0497 0.0500 0.0499 0.0506 0.0504 

CoVaRo.99'o.99(l"|X) 0.0103 0.0099 0.0104 0.0105 0.0103 

CoVaRo.95|o.99(>"|^) 0.0100 0.0099 0.0100 0.0102 0.0101 

CoVaRo.99!o.95(>"|^) 0.0501 0.0493 0.0499 0.0508 0.0507 

Table 2: Violation rates in the bivariate t(3) case. Monte Carlo backtesting 
with n = W and a, /3 e {0.95,0.99}. 




Figure 8: Bivariate t{3) samples (size n = 2000) and the joint excess regions 
in the backtesting experiment for a = P = 0.95. 

Backtesting and violation rates 

The backtesting study was implemented analogously to the bivariate Gaus- 
sian example. The results are shown in Table 2, and they go in line with 
those from the Gaussian case. While CoVaR - again, by construction - has 
a violation rate close to 1 — /3, the violation rates of CoVaR^ are significantly 
higher and increase in p. Going up to 36% for p = 0.9, the violation rates 
for CoVaR are even higher than in the Gaussian model. 

The corresponding sample plots with lines marking VaRQ,(X), CoVaR^(y |X), 
and CoVaRa{Y\X) are shown in Figure 8. Similarly to Figure 3, these graph- 
ics demonstrate how increasing dependence parameter p changes the shape of 
the corresponding sample clouds and increases the numbers of joint excesses. 
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ACoVaR as a function of p A""'''CoVaR as a function of p 




P P 

Figure 9: ACoVaR^ and A™'''^CoVaR^ as functions of p in the bivariate t(3) 
model. 

Risk contribution measures ACoVaR" and A^'^'^CoVaR" 

The comparison of ACoVaR" and A^'^'^CoVaR" is shown in Figure 9. The 
graphics demonstrate clearly how these CoVaR" based risk contribution mea- 
sures inherit the inconsistency of CoVaR". Both ACoVaR" and A'^'^'^CoVaR" 
fail to be increasing with respect to the dependence parameter p, and the 
shapes of the corresponding curves are similar to those of CoVaR" in Fig- 
ure 6. Although A^'^'^CoVaR" is slightly better behaved than ACoVaR", it 
is still strongly inconsistent with respect to p. In particular, this example 
demonstrates that the monotonicity of A^'^'^CoVaR" with respect to p in the 
Gaussian case is a special property of the bivariate Gaussian model, so that 
the advantage of A^'^'^CoVaR" over ACoVaR is quite limited. 

Extension from CoVaR to CoES 

The comparison of CoES vs. CoES" is shown in Figure 10. The monotonicity 
or non- monotonicity in p is again inherited from Co VaR or CoVaR" . See also 
Corollary 3.10(a). 

4.3 Gumbel copula with t margins 

The last model we consider here is obtained by endowing a bivariate Gumbel 
copula (cf. (10)) with t margins. Thus it has the same heavy-tailed margins as 
the previous example, but a different dependence structure. An illustration 
of the sample clouds generated from this distribution is given in Figure 13. 

On the qualitative level, all comparison results obtained in this case are 
similar to the bivariate t model, so that a brief overview is fully sufficient: 



23 



CoES as function of p 



CoES as function of p 
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Figure 10: CoES=(r|X) and CoES^(r|X) in the bivariate t(3) model as 
functions of p. 



• Corollary 3.8 guarantees that CoVaRa^/j is increasing with respect to 
the dependence parameter 6, whereas CoVaR^^^ fails to be increasing 
if dependence is at its largest (see Figure 11 for the case a = (3). The 
strongest decay of CoVaR" takes place for 6 G (1.5, 2) and slows down 
for 9 > 2. On the other hand, CoVaRa is almost constant for 6 > 2. It 
seems that for 6* > 2 the joint distribution of large values of {X, Y) is 
almost comonotonic, so that there is no much change after 6 exceeds 2. 

• The ratios CoVaRa(F|X)/ VaRa(l^) are ordered correctly with respect 
to 9, whereas the ratios CoVaR^(y|X)/ VaRQ,(y) are not (see Fig- 
ure 12). 

• The violation rates for CoVaR^^^ in a simulated backtesting study are 
significantly larger than 1 — /3, going up to 40% for a = /3 = 0.95 
and 6 = 3 (cf. Table 3 and Figure 13). This is even more than in the 
bivariate t case. 

• Both ACoVaR" and A^'^'^CoVaR" fail to be increasing in 9 (Figure 14). 

• Again, CoES is increasing in 6 and CoES^ is not; see Corollary 3.10(a) 
and Figure 15. 
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CoVaR as function of B 



CoVaR as a function of 6 
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Figure 11: Gumbel copula t(3) margins: CoVaR^(F|X) and CoVaRQ(F|X) 
as functions of Q. 



Ratio CoVaR /VaR as a function of a Ratio CoVaRA/aR as a function of a 




Figure 12: Gumbel copula with t(3) margins: Ordering of the ratios 
CoVaR=(y|X)/ VaR„(y) and CoVaR„(F|X)/ VaR„(r) for different a. 
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Gumbel/t(3) for B = 1.1 



VaRo.95 
CoVaR; 



Gumbel/t(3) for e = 1 .4 



Gumbei;t(3) for 6 = 3 




Figure 13: Gumbel copula with t{3) margins: simulated samples (size n 
2000) and the joint excess regions in the backtesting experiment for a = f3 
0.95. 
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A™°CoVaR as a function of 6 
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Figure 14: Gumbel copula with t{3) margins: ACoVaR^ and A^'^'^CoVaR^ 
as functions of 9. 
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e = 1 


6 = 1.1 


6 = 1.2 


e = 1.5 


6 = 2 


^ = 3 


CoVaRo=95,0.95(>^l^) 


0.0498 


0.0982 


0.1282 


0.1911 


0.2771 


0.4090 


CoVaRo=99,o.99(^l^) 


0.0101 


0.0346 


0.0461 


0.0752 


0.1321 


0.2423 


CoVaRo=95,0.99(>^l^) 


0.0098 


0.0309 


0.0434 


0.0754 


0.1319 


0.2450 


CoVaRo=99,o.95(^l^) 


0.0500 


0.1050 


0.1335 


0.1916 


0.2745 


0.4043 


CoVaRo.95,0.95(V'|X) 


0.0498 


0.0494 


0.0503 


0.0498 


0.0501 


0.0502 


CoVaRo.99,0.99(>"|^) 


0.0101 


0.0099 


0.0101 


0.0102 


0.0100 


0.0097 


CoVaRo.95,0.99(>^|^) 


0.0098 


0.0099 


0.0100 


0.0099 


0.0100 


0.0098 


CoVaRo.99,0.95(>^|X) 


0.0500 


0.0497 


0.0499 


0.0492 


0.0503 


0.0492 



Table 3: Violation rates for the Gumbel copula with t(3) margins: Monte 
Carlo backtesting with n = W and a, (3 e {0.95,0.99}. 



CoES as function of 




Figure 15: Gumbel copula with t{3) margins: CoES^{Y\X) and 
CoESa{Y\X) as functions of 9. 

5 Conclusions 

The present paper demonstrates that the alternative definition of Conditional 
Value-at-Risk proposed in [14, 20] (here CoVaR) gives a much more consistent 
response to dependence than the original definition used in [2, 3, 4] (here 
CoVaR=). 

The general results in Section 3 show that the monotonicity of CoVaRa,(3{Y\X) 
with respect to dependence parameters is related to the concordance ordering 
of bivariate distributions or copulas. This gives the notion of CoVaR based 
on the stress scenario {X > VaRQ,(X)} a solid mathematical fundament. On 
the other hand, comparative studies in Section 4 show that conditioning on 
{X = VaRQ,(X)} makes CoVaR" and its derivatives unable to detect sys- 
temic risk where it is most pronounced. Related counterexamples include 
several popular models, in particular the very basic bivariate normal case. 
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Based on these results, we claim that, if Conditional Value-at-Risk of 
an institution (or system) Y related to a stress scenario for another institu- 
tion X should enter financial regulation, then it should use conditioning on 
{X > VaRa(X)}. This kind of stress scenario has a much more meaningful 
practical interpretation than the highly selective and over-optimistic scenario 
{X — YaRa{X)}. Conditioning on {X > YaRa{X)} also makes CoVaR more 
similar to the systemic risk measures proposed in [15, 1, 29, 16]. 

The question how to define risk contribution measures based on stress 
events to the financial system is currently open. Besides CoVaR, CoES with 
proper conditioning may also be an option. The advantage of CoES over 
CoVaR is its coherency. In the case VaR vs. ES, this point has gained new 
interest from the regulators [7, 13]. 

In some sense, CoVaR" repeats two times the design error that is re- 
sponsible for the non-coherency of VaR. In the first step, it follows the VaR 
paradigm and thus favours a single conditional quantile of Y over an average 
of such quantiles. In the second step, it favours the most benign outcome 
of X in a state of stress over considering the full range of possible values 
in this case. Financial regulation based on CoVaR" has a strong potential 
to introduce additional instability, to set wrong incentives, and to create 
opportunities for regulatory arbitrage. 

Another argument supporting CoES is that it is particularly suitable for 
stress testing. In a system with several factors Xi, . . . , X^, the numbers 
CoESai,i3{Y\Xi) describe the influence of the different Xi on Y. Assigning 
relative weights Wj to the scenarios X^ > VaRo,. {Xi) and taking the weighted 
sum 



one always obtains a sub-additive risk measure. If the weights Wi sum up to 
1, the resulting risk measure is coherent in the sense of [5]. The choice of the 
weights Wi or of the confidence levels may change over time, incorporating 
the newest information about the health of the institutions Xi, . . . ^Xd- 

To make the weighted risk measure (18) even more meaningful, one could 
modify it by implementing not only the single risk factor excesses Xi > 
VaRa,.(Xj), but also the joint ones. Consistent choice of the corresponding 
weights can be derived by methods presented in [23]. A detailed discussion 
of this goes beyond the scope of the present paper and would also require 
additional mathematical research. 

Motivated by the recent financial crisis and the following discussions on 
appropriate reforms in financial regulation, systemic risk measurement has 



d 





1=1 
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become a vivid topic in economics and econometrics. Our results show that 
some important contributions are also to be made in related mathematical 
fields, including probability and statistics. In particular, the dependence 
consistency or, say, dependence coherency of systemic risk indicators is a 
novel problem area that needs further study. The present paper provides first 
examples and counter-examples for compatibility of systemic risk indicators 
with the concordance order. The questions for general characterizations or 
representations of functional with this property are currently open. 

In addition to dependence consistency, implementation of systemic risk 
measures in practice obviously needs estimation methods. The estimation of 
CoVaR in GARCH models is discussed in [14]. As non-parametric estimation 
of rare events would needs a lot of data, methods from Extreme Value Theory 
may be used to extrapolate the rear events from a larger number of data 
points. A similar approach for conditional default probabilities is pursued 
in [29]. 
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